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Abstract 



We analyze the expectation value of observables in a scalar theory on the fuzzy two sphere, 
represented as a generalized hermitian matrix model. We calculate explicitly the form of 
the expectation values in the large- iV limit and demonstrate that, for any single kind of 
field (matrix), the distribution of its eigenvalues is still a Wigner semicircle but with a 
renormalized radius. For observables involving more than one type of matrix we obtain a 
new distribution corresponding to correlated Wigner semicircles. 



1 Introduction 



Matrix models have been of interest for a long time going back to Wigner's work denning 
the classic Gaussian ensembles for matrices pQ . The motivation for this was to understand 
the distribution of energy levels for heavy nuclei, for which the Hamiltonian has so many 
contributing interactions that it my be taken to be a random NxN matrix. The distribution 
of eigenvalues in the large iV limit could then be compared to the distribution of energy 
levels for a large collection of different nuclei. Assuming a normal distribution for the 
elements of the matrix, the classic semicircle law of Wigner emerged in this context. Since 
then matrix models have emerged in many other contexts in physics: in modeling Riemann 
surfaces with a view to applications in string theory [S], in integrable systems, such as the 
Calogero model, as well as other condensed matter systems [3], in possibly elucidating the 
concept of chaos in quantum systems 0] , etc. 

Since matrices are a simple example of noncommuting variables, Wigner's matrix en- 
sembles also turn up in the noncommutative probability theory of Voiculescu [5j. In fact, 
the semicircle law plays a very important role in free (uncorrelated noncommutative) prob- 
ability theory. Essentially it is to free probability theory what the normal distribution is to 
the probability theory of commuting variables, with a corresponding central limit theorem 

An obvious question that arises is whether there are other matrix ensembles of impor- 
tance in physics and mathematics which are also naturally defined. Fuzzy spaces are almost 
a self-evident answer to this question. Noncommutative (fuzzy) spaces, and field theories on 
such spaces, have been an important topic of research for a long time now [6]. Such spaces 
can arise as brane solutions in certain contexts in string theory and in the matrix version of 
M-theory [7] . Gauge theories on such spaces are interesting since they can describe fluctua- 
tions of the brane solutions and unify in a natural way gauge and spatial degrees of freedom. 
This has generated interest more generally on field theories on fuzzy and noncommutative 
spaces. 

Fuzzy spaces are noncommutative spaces that can be described by finite dimensional 
matrices, and, by now, there are many examples of such spaces. When the dimension of 
the matrices becomes large, these spaces tend to corresponding commutative manifolds in 
terms of their geometry and the algebra of functions on such spaces. From the physics point 
of view fuzzy spaces are important for many reasons: 1) They provide a regularization that 
can preserve various symmetries, even avoiding the fermion doubling problem as compared 
to the standard lattice regulator. 2) They do play a role as effective descriptions of cer- 
tain condensed matter systems such as the quantum Hall effect. 3) Being a finite- mode 
approximation to fields preserving isometries, they have implications for quantum gravity. 
4) Two dimensional YM theory naturally reduces to a (unitary) matrix model with time as 
a continuous parameter. 

There are, of course, many diverse issues here, but if we take the simple case of a field 
theory on a fuzzy space, with the fields being NxN matrices, the Euclidean functional 
integral is a matrix ensemble. The Gaussian ensembles of Wigner correspond to mass terms 
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for the fields. The action also has kinetic energy terms, given, for a scalar field, by the matrix 
Laplacian. The Laplacian essentially defines the geometry of the continuum manifold as 
we take the large N limit and hence it is the key geometrical ingredient. This clearly gives 
a natural class of matrix ensembles. So, unlike the Wigner distribution, we must seek a 
large ./V limit which takes account of the contribution of the Laplacian. We also have new 
observables in this case and the Laplacian will give some degree of correlation which is 
reflective of the emergent geometry at large N. 

Not surprisingly, there have been attempts to understand the role of and generalize the 
Wigner distribution for fuzzy spaces The authors of |Sj, in particular, start with the 

mass term as the leading term of the action and integrate out the angular modes which 
appear in the kinetic term, treating this term as a perturbation. There is, however, no 
guarantee of continuously connecting this to the zero mass limit because of possible phase 
transitions as a function of N and the lack of summability of the perturbation series. Hence 
an approach treating all terms of the action on an equal footing with the possibility of 
different scaling limits is desirable. 

The formulation of scalar field theory on the fuzzy sphere by Steinacker, on the other 
hand, is closer to the spirit of our work. Using momentum space integration, rather than 
matrix techniques, the Wigner distribution for the eigenvalues of the matrix is recovered 
[9]. The main difference with the present work is our focus on observables containing also 
derivatives of the field (Laplacians and their powers). These manifest the dependence of the 
results on the angular variables of the matrix and introduce nontrivial correlations. The 
calculation of mixed expectation values, with both powers of the field and its derivatives, 
and the resulting correlated Wigner distribution are our main results. 



2 Ensembles and distribution functions 

2.1 The ensembles 

We start with the simplest case of a real scalar field on a fuzzy two-sphere. The Cartesian 
coordinates are N x N matrices: 

X a = —L a , [L a ,Lp] = ie a p 7 Lry , y^X%= (l - j r 2 (1) 

OL 

where L a , a = 1,2,3 are SU(2) generators (angular momentum matrices) in the N- 
dimensional representation, r represents the radius of the sphere and 9 = 2r 2 /N the noncom- 
mutativity parameter. Fields become general N x N matrices M. Derivatives (rotations) 
and the corresponding Laplacian are L-commutators 

C a M = -i[L a , M] , AM = -i J2l L «> ^ M H ( 2 ) 

a 

while integration over the sphere becomes a matrix trace 

C 4:7rr^ 

/ d 2 x<& = tiM (3) 

Js* N 
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A real free scalar field on the two-sphere is represented by a hermitian N x N matrix M. 
Upon rescaling M to normalize the kinetic term, the Euclidean action is given by 

S = ~Tr([L a ,M][L a ,M}) + ^Tr(M 2 ) 

= ^Tr(M [L a , [L a , M]]) + ^Tr(M 2 ) (4) 

We shall be interested in the large- iV limit of Q. Classically, for 'smooth' configurations, 
M is replaced by its symbol in this limit, which is now a real scalar field 4> on the two-sphere. 
The adjoint action of L a on M, i.e. [L a ,M], becomes the gradient of cj) and the action, 
upon appropriate scaling, becomes the free field action on the two-sphere, 

S = \f d^(S 2 ) [(V0) 2 + /iV] (5) 

Quantum mechanically, however, expectation values of observables depend on the noncom- 
mutativity parameter, which acts as a regulator of infinities in the continuum, and the 
large- N limit remains nontrivial. 

For the usual Gaussian matrix ensembles for which we just have the mass term in 
Q, we have the symmetry M — > U^MU for a unitary matrix U, which allows us to 
diagonalize M as M = W M^iagU . Only the eigenvalues of M are physically meaningful 
and their distribution is given in the large ./V" limit by the Wigner distribution function. 
The observables of interest are of the form TrM fc . In the case of Q, it is still possible 
to scale variables and the parameter [i in such a way that the Laplacian of M, which we 
denote by B = [L a , [L a , M]], becomes irrelevant in the large iV limit and one recovers just 
the Wigner result. However, we are interested in cases where the gradient term is not 
irrelevant, particularly for possible application to field theory on fuzzy spaces. In fact, we 
will consider an action which is more general than the one given in Q, 

1 a 2 

S=-Tt(MKM) + ^-Tr(M 2 ) (6) 

where the kinetic operator K can be more general than (Lq ) 2 . 

The observables involve arbitrary products of M, B = KM, [L a ,M], etc. Rather 
than changing variables to the eigenvalues (M^iag) and the angular degrees of freedom 
(U), we will consider the moment generating functions for the observables and obtain the 
distributions functions which lead to them. For this, the form of the action itself will 
not be very important; in fact ^ is meant more as a guide or motivation. All we need 
is that correlators can be evaluated by Wick contractions in terms of products of two- 
point functions. Rather than the action, we can specify the scaling behavior of two-point 
functions. We will examine possible scalings which give a finite large N limit. 

2.2 Recursion rules and the generating function for correlators 

To facilitate the upcoming calculations, we define a basis for N x N matrices {T^}, where 
A = 1, 2, • • • ,21 + 1 for each I and / = 0, 1, • • • , N — 1. These matrices rl transform as the 
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spin-Z representation of SU(2). Explicitly, a possible choice for the elements of this base is 
1, L a , L a L/s + LpL a — (2/3)L 2 8 aj 3, etc. We take T A to be normalized as 

Tr (TfT { P)=5 11 ' 5 AB (7) 

Further, [L a ,[L a ,T^]] = 1(1 + The matrix M may be expanded in this basis as 

M = ^ A c l A Ti . The action is diagonal in terms of this basis, with the expectation value 
(c l A c l B ) = SabS 11 ' G{1). For the action Q, G(l) = [fj 2 + 1(1 + l)]" 1 , but it can be taken to 
be a more general function of I depending on the form of the kinetic term operator K in 
Using the fact that the matrices T A transform as the spin4 representation of SU(2), 
we see that 

,(/) (J )\ _ 21 + 1 



A 

This follows easily from carrying out an SU(2) transformation on the states corresponding 
to the matrix labels i, j. 

First consider the propagator for M, or the two-point function (MijM^i). As a result 
of relation ([8]) , the two-point function (or propagator) for a matrix product of two matrices 
takes the form 

1 N ^ 

( (M M)y> = ^ Yl & + l ^ 5 v = f ( 9 ) 

1=0 

Similarly, we can define the propagators 

((5%) =9^, ((MB) t ,) = ((BM)ij) = h5 lJ (10) 
For B = [L a , [L a , M]], we can explicitly obtain 

N-l N-l 
9=nT,( 21 + V l2 ( l + l ) 2 °^ h =N^ {2l + m + 1)G{1) 

1=0 1=0 

Note that, for positive G(l), an application of the Schwartz inequality to the above sums 
gives the relation 

fg > h 2 (12) 

An important result in the large limit is that, for any matrix product, the expectation 
value is proportional to the identity matrix. To see this, consider (Aij) where A is a matrix 
of the form 

A = M mi B hl M m2 B b2 ■ ■ ■ (13) 

In evaluating the expectation value of A^ by Wick contractions, the leading term in the 
large N limit will involve only planar contractions; i.e., diagrams where the the propagators 
"cross" are subleading in 1/N. (This can be shown by examining the scaling of the matrix 
elements for J"! that appear in such contractions and the corresponding factors of N arising 
from index summations.) As a result, in the leading term, there will be at least one case 
where two adjacent matrices are contracted. Denoting Q1Q2 the two adjacent matrices, 
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and L and 1Z the matrices to the left and right of Q1Q2, and using ( 9|10 ), this contraction 
will give us a term of the form 

(Aij) « (C ia TZ bj ) ((Q1Q2U) ~ (Wy) 9 (14) 

where Q\, Q2 stand for either M or 5 and q is the appropriate function among /, g, h. 
This step reduces the correlator to one with lower number of matrices M or B, with the 
two nearby matrices that were Wick contracted deleted. Iterating, we see that 



(A 



a5i 



where a = (trA) /N is a scalar quantity. 



(15) 




m b 



m-1 b-a-1 



Figure 1: Illustration of equation ( 16 ). Circles represent expectation values of matrices. The 
thin line represents M's, the solid line B's and the dashed line represents the contraction. 
Planarity ensures that there are going to be no contractions between the two parts. 

We now turn to the correlator (Tv(M m B b )}. In evaluating this by Wick contractions, 
first consider the contraction of an M matrix next to the series of -B's with one of the B's. 
The result may be written as (see Figure 1) 

1 



<(M 



m—l\ 



Mjk {B a )kr B rs (B 



6-a-l\ 



h „(TrB< 



(Tr(M m - 1 B b - a - 1 )) 



(16) 



In arriving at this, we have used two results: Since nonplanar contractions are suppressed 
at large N, the contractions of the series of -B's in B a have to be within themselves, and 
hence the series can be replaced by the expectation value. Further, since the expectation 



value is proportional to the identity from (15), the remaining set of matrices fall into a 



{{M% M jk (M a ) kr M rs (M 



m-2-a-P 



matrix product. Finally, the contracted matrices, indicated by the overdots, give a factor of 
h. Similarly, considering the contraction of the chosen M matrix with another M matrix, 
we get 

B b )si) = f ^(TvM a ) (Tr(M m - 2 - a B b )) (17) 

Combining (16) and (17) and allowing for all possible values of a, we get the recursion rule 

m-2 

(Tr(M m B b )) = fjrJ2 ( TrM °) (Tr{M m - 2 - a B b )) 

1 b ~ l 

+ h — ^(Tr(M m - 1 B b - a - l ))(TrB a ) (18) 



a=0 
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Since one M is chosen to be contracted with all other matrices, this relation applies when 
m > 0. 

In a similar way, we can single out a B matrix adjacent to the series of M matrices and 
consider its contractions. This leads to the recursion rule, for b > 0, 

1 b ~ 2 

(Tr(M m B b )) = g — Y^ (TrB a ) (Tr(M m B b ~ 2 ~ a )) 

1 m— 1 

+ h — Y (Tr(M m_a_1 B b - 1 )) (TrM a ) 



(19) 



o=0 



These recursion rules may also be viewed as the Schwinger-Dyson equations for expectation 
values calculated via the functional integral 



(O) = / [dM] e~ s ^ O 



Equation ( 19 ) can be obtained, for example, by considering the identity 

d 



[dM] 



dMji 



(M m_1 B b ~ l )ji e 







(20) 



(21) 



We will not consider the simplifications of such matrix integrals, but rather proceed to the 
direct solution of the recursion rules. For this, we now define the normalized correlator 



W m>b = i / Tr 



M 



B 



(22) 



The two recursion rules become 

m-2 



6-1 



4 W m , b = Y W a , W m -2- a> b + 7 Y W ' 



rn-l,b-a-lWo a 



a=0 
b-2 



a=0 
m—1 



4 W mjh = Y W , a W m)b -a-2 + 7 Y W m-a-l,b-lW ai0 



a=0 



a=0 



where 
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Define the generating function 



<t>(t,s)= Y W m , b t m s b 



(23) 
(24) 
(25) 

(26) 



m,6=0 



which is normalized as 4>(0, 0) = Wq^q = 1. The recursion rules now become 

4[(f>(t,s)-(f>(0,s)} = t 2 (f)(t,s)(t>(t,0) + jts</)(t,s)<t)(0,s) 
4 [</>(t, s) - (j)(t, 0)] = s 2 0(t, s) <M0, s) + -yts<t>(t,s)<(>(t,0) 



(27) 
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Solving these equations for (f)(t, s) and equating the two expressions, we find 

4(j)(t, 0) - t 2 <j)(t, 0) 2 = 40(0, s) - s 2 0(O, s) 2 (28) 

Thus each expression must be a constant, which should be 4 from 0(0,0) = 1. We then 
solve ([28| and use it in (|27|) to get 



<l>(t,s) 



l + x/r^t 2 ) l + v 7 ! 



7^s 



(29) 



Also, as a consequence of (12), we have for the constant 7: 

- 1 < 7 < 1 



(30) 



The value of 7 is crucial. For 7 = 0, in particular, the above generating function is the 
product of the generating functions of two independent Wigner distributions, while for 
nonzero 7 we have correlations. 



2.3 Scaling and special cases 



To explore the possible values of 7, it is useful to consider the scaling properties of the 
functions /, g, h. Going back to their definition in (9pT), and taking G(l) to go like l a for 
large values of I, we see that 



(31) 



N ^ N 

The summations can be approximated by integrations when the exponent is larger than 



— 1; otherwise they are dominated by small values of I. From the behavior in (31), we see 
that there are four possible cases. 



1. a > -2 : 

In this case / ~ iV Q+1 , g ~ iV Q+5 , h ~ N a+3 and so, 7 ~ 1 (order N°). 

2. -4 < a < -2 : 

In this case / ~ iV -1 , g ~ N a+5 , h ~ N a+3 and so, 7 ~ iV2( a+2 ) < 1. 

3. -6 < a < -4 : 

This leads to / ~ A^" 1 , g ~ N a+5 , h ~ JV _1 , and hence, 7 ~ jV~2( a + 6 ) < 1. 

4. q < -6 : 

In this case, / ~ N , g ~ N , h ~ A^ _1 and hence, 7 ~ 1. 



In cases 2 and 3, the mixed term will be irrelevant (7 —> 0) and 4>(t, s) becomes the product 
of the moment generating functions for two independent Wigner distributions. In the other 
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two cases, 7 ~ 1, and the distribution will be different from the Wigner semicircle law. The 
case of a = —2 is special. In this case, 7 vanishes as (lniV) -1 . The distribution function 
will have to be evaluated by taking a proper limit as 7 tends to this value. Notice that 
a = —2 is the case for the kinetic term being the Laplacian. A similar statement applies to 
a = —6. 



2.4 The distribution function 



In oder to obtain the distribution function from the moments, we start with (29) and expand 
4>(t, s) in 7 as 



n=0 



= £ Z t n 4>{t) n+1 s n 4>{s) n 



+1 



where 



(«) = ^(«,0) = ^(0,«) 



l-v 7 !^? 



1 + Vl -u 2 



So to each order in 7 the distribution factorizes as 

00 n 

pO, y) = Yl (z) Pn(v) 



n=0 



(32) 
(33) 

(34) 



To find the n-th order distribution p n (x) corresponding to the generating function t n (p(t) n+1 
we work as follows: Define [/]+ the non-negative power (non-singular) part of a function 



[ h a- 2 t 2 + a-it 1 + a + a\t + a 2 t 2 + •••], = a + ait + a 2 t 2 + 



(35) 



Then, if the generating function f(t) corresponds to the distribution w(x), the generating 
function 

7W" 



t n 



, n > 



x n w(x) 



(36) 



corresponds to the distribution x n w(x), as can easily be shown. (Note that this is not true 
for n < 0.) 

Now, it is easily verified by direct substitution that the function 4>(t) obeys the identity 



4>{t) 



n+l 



1 



1 



n+l / , \ n+l 

1 + Vi-t 2 ) -(1-Vi-t 2 



4>(t) 



1 + Vi-t 2 ) -(1-Vi-t 2 



VT^ 2 

Upon expanding the binomials and multiplying with t n we obtain 



(37) 



2\A: 



fc=0 



(38) 
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We observe that the coefficient of <fi(t) contains only negative or zero powers of t, while 
the second term contains only negative powers. The left hand side, however, is obviously 
nonsingular. Therefore, the second term cancels the singular part of the first term but does 



not contribute to the nonsingular part. Using the result (36) we obtain 



[n/2] 



Pn{x) = ^ 



k=0 



n + 1 
2k + 1 



2 n x n - 2k {x 2 - l) k p{x) 



(39) 



where p(x) is the Wigner distribution, 



p(x) = -\pY 

IT 



(40) 



We can now consider the two-dimensional distribution p(x, y) as defined by 

W mb = / dxdy p(x,y) x rn y b 



(41) 



The distribution p(x,y) can be calculated by substituting for p n from (39) in (34) and 



carrying out the summation. This task is facilitated by noticing that (39) can be obtained 
from the first term in (37) by substituting t — ' 



p(x,y) = p{x)p{y) 



x . After some algebra we find 
1-7 2 



(1 - 7 2 ) 2 - 47(1 + j 2 )xy + 47 2 (x 2 + y 2 ) 



(42) 



This is symmetric in x and y. The integral of p(x, y) over y is nontrivial; it can be calculated 
and gives p(x) as it should, giving a check on the expression for p(x,y). 

Although the quadratic form in x, y in the denominator is not positive definite, the 
above distribution is positive for 7 2 < 1. Indeed, for 7 > the negative eigenvalue of the 
quadratic form corresponds to the eigenvector x = y. Putting x = y = 1, their maximal 
value given the Wigner distribution prefactor p(x)p(y), the denominator becomes (1 — 7) 4 . 
Similarly, for 7 < the negative eigenvalue corresponds to x = —y and choosing x = —y = 1 
the denominator becomes (1 + 7) 4 . For 7 2 = 1 the distribution appears to be singular, but 
taking the limit we see that it becomes p{x)5{x — y) for positive 7 and p{x)5{x + y) for 
negative 7. 

For 7 2 > 1, however, the distribution becomes negative, signaling the nonexistence of 
a probability interpretation in that case. We can understand this by noticing that the 
correlation of the variables x and y is calculated as 



(xy) 



2d t d s 



7 



(43) 



t=s=0 



Since correlations of stochastic variables have to be between —1 and 1, this is also the 
allowed range of 7. For 7 = ±1 the two variables are fully (anti)correlated and thus equal 
to (minus) each other, justifying the delta-functions. 
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The relation (30) ensures that the distribution is always positive. For instance, since 7 
is the ratio h/yffg, in the regular large- iV scaling case corresponding to G(£) ~ N a , a > —2 
(case 1 of subsection 2.3), it becomes 

= V(« + 2)(a + 6) 

We see that 7 is always between and 1, and it tends to 1 in the limit a — > 00. In this limit, 
it is reasonable that the eigenvalues of M and B are correlated since they are dominated 
by the largest angular momentum sector. 

In conclusion, we obtain a correlated distribution for the eigenvalues of M and B, while 
their marginal distributions remain Wigner semicircles of radii 2y/J for M and 2^/g for B, 



as seen from the scaling factors in (22). 



3 The massless case 



The action describing a massless scalar on the fuzzy sphere, \x = in the action Q, 
merits special attention. This is the critical case a = —2 in the scaling of G(l) at large I, 
corresponding to a very weak vanishing of 7 ~ (In N)~ 1 / 2 at the large N limit. Further, for 
\i = the constant mode of the field (the trace of M) drops out of the action and must be 
eliminated from the calculation. 

In fact, taking the limit [i — > in this case is somewhat nontrivial. The propagator is 
G(l) = [fi 2 + 1(1 + and for fi ~ 1/lniV or smaller we obtain in the large N limit 

/ = i( M - 1 + 21niV) , g = \N\ h = N (45) 

which leads to a Wigner radius of 2yff = 2\J ([i~ l + 2 In N) /N for the eigenvalue distribu- 
tion of M. This not only diverges as /i goes to zero, but is also misleading: for such low 
values of /U the planarity property of matrix expectation values fails, since the trace part 
of M contributes to the same order (or higher) than the traceless part and arises in all 
diagrams (planar and nonplanar). Such contributions give rise to a Gaussian, rather than 
Wigner, distribution. 

To understand this better, we decompose M in its trace part cq and its traceless (I 7^ 0) 
part M 

M = 4L+M = ^L+y c l A T^ (46) 

The trace and traceless parts decouple. The eigenvalues of M have a Wigner distribution at 
large N with radius 2^/2 In N/N, arising from / with I = dropped, while the single mode 
of the trace part contributes a shift distributed as a Gaussian with spread 1 / \fWJjL. Since the 
two distributions are independent, the total eigenvalue distribution will be given by their 
convolution. For \i 3> 1/ In N the spread of the Gaussian is much smaller than the Wigner 
radius and the convolution essentially gives back the Wigner. For fi <g 1/ In JV, on the other 
hand, the Gaussian dominates. For n ~ 1/lnA^ we get an intermediate distribution. 
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In order to obtain the massless result, therefore, we have to omit the trace part of the 
matrix (which has vanishing action) and substitute M for M. The result is again a Wigner 
semicircle of radius y/8 In N/N for M, correlated weakly (7 ~ (IniV) -1 / 2 ) with a Wigner 
semicircle for B of radius = V2N 3 . To first order in 7, odd-odd expectation values of 
the form W^2fc+i,2o+l a l so survive and are given by (7/4)^1 (x) pi(y) as in (34). 



4 Conclusions and discussion 

We recovered a correlated distribution that reproduces the expectation values of (ordered) 
matrix products of the basic matrix variable M and its Laplacian B in the large- N limit. 
This is our main result. 

The distribution of eigenvalues of M remains a Wigner semicircle, but with a renor- 



malized radius R = 2\J~j as is evident from (22) and (42). Similarly, the distribution of 
eigenvalues of B is also a Wigner semicircle with radius R' = 2^fg. In this sense, the Wigner 
semicircle is very robust, arising basically from the planarity property of expectation values 
of matrix observables in the large- N limit. This result is at odds with the result of [H] which 
gives a polynomially deformed Wigner distribution in the presence of the kinetic term in the 
action. The source of the discrepancy could be the fact that that the kinetic term becomes 
dominant in the large- ./V limit, while it was treated perturbatively in [8]. 



The new element in our result is the correlation 7. The distribution (42) is a two- 
dimensional correlated Wigner distribution (whose one-dimensional marginals are the stan- 
dard Wigner semicircle). Indeed, the expectation values of mixed products of M and B do 
depend on the angular variables of M and these introduce the nontrivial (and non-unit) 
correlations. This is independent of the exact form of the action, as long as the trace sector 
does not dominate and 7 remains nonzero in the large- iV limit. For instance, for the stan- 
dard matrix model (without the kinetic term), that is, the case a = in the scaling of G(l), 
we have 7 = V3/2. For the case with the kinetic term, a = —2, the correlation vanishes 
weakly (logarithmically) in the large- N limit. 

There are obviously many issues remaining to be analyzed; namely the expectation value 
of more general products of matrix variables, the question of the dependence of the results 
on the ordering of matrices and the relevance of the results to string, membrane or gravity 
models. These will be examined in forthcoming publications. 

Further, we point out that our distribution represents an explicit expression of correlated 
free variables in Voiculescu's noncommutative probability theory. Large- N matrix models 
with general propagator G(l), and their generalized matrix observables, are a natural arena 
where these mathematical notions are realized and make their physical significance clear. 

The direct physical relevance of our results lies mainly in noncommutative theory and 
its properties, contrasted to those of a regularized standard (commutative) field theory. 
We see that the noncommutative theory presents qualitatively different features, mainly 
related to planarity, which is itself an expression of UV/IR mixing. The generic emergence 
of Wigner, rather than gaussian, distributions, and their correlated generalizations, is the 
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key signature of this effect. Other possible physical applications of our results range from 
quantum Hall situations to quantum gravity. Matrix models for the quantum Hall effect 
have been proposed and used to some advantage [10], and the inclusion of a kinetic term 
to such models adds an element of compressibility that could be probed by an approch 
similar to the one in this work. Finally, quantum gravity remains the main motivation 
behind this work and the possibility to probe Planck-scale effects, such as the formation 
and evaporation of microscopic black holes, without the benefit (or burden) of string theory 
is an exciting prospect. 
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